BIFURCATION DIAGRAM AND STABILITY FOR 
A ONE-PARAMETER FAMILY OF PLANAR VECTOR FIELDS 

J. D. GARCIA-SALDANA, A. GASULL, AND H. GIACOMINI 

Abstract. We consider the 1-parameter family of planar quintic systems, x = y s —x i , 
y = —x + my , introduced by A. Bacciotti in 1985. It is known that it has at most 
one limit cycle and that it can exist only when the parameter m is in (0.36,0.6). 
In this paper, using the Bendixon-Dulac theorem, we give a new unified proof of all 
the previous results, we shrink this to (0.547,0.6), and we prove the hyperbolicity of 

fvj \ the limit cycle. We also consider the question of the existence of polycycles. The 

main interest and difficulty for studying this family is that it is not a semi-complete 

f>) family of rotated vector fields. When the system has a limit cycle, we also determine 

explicit lower bounds of the basin of attraction of the origin. Finally we answer an 

open question about the change of stability of the origin for an extension of the above 

00 ' systems. 
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1. Introduction and main results 



Q 

A. Bacciotti, during a conference about the stability of analytic dynamical systems, 
held in Florence in 1985, proposed to study the stability of the origin of the following 
quintic system 

f x = y 3 -X 3 , ,_ 

J> , | y = —x + my 5 , m G K. 



Two years later, a quite complete study of ([I]) was done by Galeotti and Gori in |10j . 
They prove that, when m G (— oo,0.36] U [0.6, oo), system ([1]) has no limit cycles and, 
otherwise, it has at most one. Their proofs are mainly based on the study of the stability 
t^J" \ of the limit cycles, controlled by the sign of its characteristic exponent, together with a 

transformation of the system using a special type of adapted polar coordinates. Their 
proof of the uniqueness of the limit cycle does not provide its hyperbolicity. 

In this paper we refine the above results. To guess which is the actual bifurcation 
diagram we did first a numerical study, obtaining the following: it seems that there 
/\ , exists a value m* > 0, such that: 

(i) System ([1]) has no limit cycles if m G (— oo, m*] U [0.6, +oo). Moreover, for m = m* 
it has a heteroclinic polycycle formed by the separatrices of the two saddle points 
located at (±m — 1 ' , dbm — 1 ' ). 
(ii) For m G (m*,0.6) the system has exactly one unstable limit cycle, 
(iii) The value m* is approximately 0.560115. 

Recall that a polycycle is a simple closed curve formed by several solutions of the system 
and admitting a Poincare return map. The first two items coincide with the ones 
described in [10]. In that paper it is claimed that m* is between 0.58 and 0.59, but 
our computations give a different result, which we believe that is the right one. 
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The first aim of this work is to obtain analytic results that confirm, as much as 
possible, the above description. To clarify the phase portraits of the system, we will 
draw them on the Poincare disc, see [3| I24j . 

For m < 0, system ([1]) has no periodic orbits because x 2 /2 + y 4 /4 is a global Lyapunov 
function. Therefore, the origin is a global attractor. In particular, its phase portrait is 
trivial. Therefore, we will concentrate on the case m > 0. In this case, the system has 
three critical points, (±m -1 ' 4 , ±m -1 ' 4 ) and (0, 0). The first couple of points are saddles 
and the third one is a monodromic nilpotent singularity. Its stability can be determined 
using the tools introduced in [21 [19], see Subsection [2] and next Theorem 11.31 We prove: 

Theorem 1.1. Consider system ([1]). 

(i) It has neither periodic orbits, nor poly cycles, when m £ (— oo, 0.547] U [0.6, oo). 
Otherwise, it has at most one periodic orbit or one polycycle and both can not 
coexist. Moreover, when the limit cycle exists, it is hyperbolic and unstable. 
(ii) For m > 0, their phase portraits on the Poincare disc, are given in Figure{Jl 
(iii) Let M. be the set of values ofm for which it has a heteroclinic polycycle. Then M. 
is finite, non-empty and it is contained in (0.547,0.6). Moreover, when m G M, 
the corresponding system has no limit cycles and its phase portrait is given by 
Figure [7] (b) . 





(a) When m G (0,0.547], or when 

m G (0.547, 0.6) and neither the 

polycycle nor the limit cycle exist. 



(b) When m G (0.547,0.6) and 
the polycycle exists. 





(c) When m G (0.547, 0.6) and 
the limit cycle exists. 



(d) For m G [0.6, oo) 



Figure 1. Phase portraits of system (pQ). 



Our simulations show that (a), (b) and (c) of Figure Q] occur when m G (0,m*), 
m = m* and m > m* , respectively, for some m* G (0.547,0.6), that numerically we 
have found to be m* ~ 0.560115. We have not been able to prove the existence of this 
special value m*, because our system is not a semi-complete family of rotated vector 
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fields (SCFRVF) and this fact hinders the obtention of the full bifurcation diagram; see 
the discussion in Subsection 13.11 and Example 17.11 This is precisely the reason for which 
we have decided to push forward in the study of system ([1]) . Our approach can be useful 
to understand other interesting polynomial systems of differential equations that have 
been considered previously; see for instance [3J [8] . 

In any case, from our result, we know the existence of finitely many values m*, 
j = l,...,k, where k > 1, satisfying 0.547 < m* < m\ < ■■■m^ < 0.6, such that 
phase portrait (b) only happens for these values. Moreover, for m G (0.547, m*), phase 
portrait (a) holds, for m G (m^.,0.6) phase portrait (c) holds, and for each one of the 
remainder k — 1 intervals, the phase portrait does not vary on each interval and is 
either (a) or (c). 

As a byproduct of our approach we can also give explicit algebraic restrictions on the 
initial conditions to ensure that the solutions starting at them tend to the origin. 

Recall that when a critical point, p G W 1 , of a differential system is an attr actor we 
can define its basin of attraction as 






{xG 



lim (/?(£, x) 
t—>+oo 



p}, 



where (p denotes the solution of the differential system such that (/?(0, x) = x. A very 
interesting question, mainly motivated by Control Theory problems, consists in obtain- 
ing testable conditions for ensuring that some initial condition is in Wl Usually these 
conditions are obtained using suitable Lyapunov functions. We prove next result using 
a different approach based on the construction of Dulac functions. 

Proposition 1.2. LetW^ be the basin of attraction of the origin of system ([1]). Consider 
V m (x, y) = go,m(y) + 9i,m{y)x + 92,m,{y)x 2 , with 

92,m{y) =^7T7^t( 3 - 10m)(3 + 35m)y 12 - ^^(75 - 125m) 2/3 (3 - 13m)y 8 



89100 
1 



6300 



1 



+ —(3 - 10m)y 6 - ^(75 - I25m) 2 / 3 y 2 + 1, 
90 lb 

9i,m(y) = 92, m (y) and 9o,m(y) = g'imiy)/ 2 ~ m y 5 92, m (y)/ 2 + Wffi^tl/) A Then, for 
m G (0.5,0.6), U m C Wq, whereU m is the bounded connected component of{(x,y) G M? : 
Vm(x,y) < 0}, that contains the origin and whose boundary is the oval ofV m (x,y) = 0, 
see Figure^ 




Figure 2. The limit cycle of system ([I]) and the set U n 
Proposition 11.21 when m = 0.57. 



introduced in 



As we will see, the proof of the above proposition is an straightforward consequence 
of Proposition 15.21 Using the same tools, it can be shown that the same result also holds 
for smaller values of m. In any case, notice that this proposition covers all the values of 
m for which the system has limit cycles. 
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Studying the stability of the origin of system ([1]) we realized that, using the same 
tools, we could solve an open question left in [TO]. Our third result studies the stability 
of the origin of the following generalization of system ([1]) : 

x = y 3 -x 2k+1 , , 2 * 

y = —x + my 2s+1 , m G E and k,s£N + . 

In [TO], the authors gave the stability of the origin when s ^ 2k and ask whether it 
is true or not that the change of stability of the origin when s = 2k is at the value 
m = {2k + 1)/(4A; + 1). We will prove that their guess was not correct for k > 1. Next 
result shows that when s = 2k, the stability changes at 

(2fc + l)!! 
m= (4/c + l)!!!!' (3) 

where, given re G N + , n!! and n!!!! are defined recurrently, as follows, 

re!! = n x (re - 2)!!, re!!!! = re x (re - 4)!!!!, 

with l!! = 1,2!! = 2 and j!!!! = j for 1 < j < 4. Notice that when k = 1, the right 
hand-side of ([3]) and (2k + 1)/(4A; + 1) coincide and give m = 3/5, which is one of the 
values appearing in Theorem 1 1.11 

Theorem 1.3. Consider system ([2]). Then: 

(i) When s < 2k, the origin is an attractor when m < and a reseller when m > 0. 
(ii) When s > 2k, the origin is always an attractor. 

(hi) When s = 2k, the origin is an attractor when m < (2k + l)!!/(4/c + 1)!!!! and a 
repeller when the reverse inequality holds. Moreover, when k = 1 and m = 3/5 
the origin is a repeller and for m < 3/5 system §Q) has at least one limit cycle 
near the origin. 

The method used to study the stability of the origin of ([2]), when s = 2k and k = 1, 
also works for deciding its stability for the cases not covered by the above theorem: 
s < 2k and m = 0; and s = 2k, k > 1 and m as in Q. Nevertheless, the computations 
are tedious and we have decided do not perform them. 

The paper is structured as follows. In Section [2] we prove Theorem 11.31 Section [3] 
collects some preliminary results. It starts with a discussion on the differences between 
being or not, a SCFRVF. Then, subsection 13.21 is devoted to study the singularities at 
infinite of system ([1]) and their phase portraits on the Poincare disc. Afterwards, we 
present some Bendixson-Dulac type results that we will use to prove non-existence or 
uniqueness of periodic orbits or polycycles. Finally, we introduce a result for controlling 
the number of roots of 1-parameter families of polynomials and we show that our system 
can be reduced to an Abel differential equation. 

In Section 0] we prove the non-existence results for m G (— oo,0.36] U [0.6, oo). Our 
proof is different of that of [TO] and it is mainly based on the use of Dulac functions. 

In Section[S]we prove the existence of at most one periodic orbit when m G (1/2, 0.6). 
Our approach also gives the hyperbolicity of the orbit and again uses a Bendixson-Dulac 
type result. This section also includes the proof of Proposition 11.21 

Section is devoted to enlarge the region where we can assure the non existence of 
periodic orbits and polycycles, proving this for m G (9/25,0.547]. The proof uses once 
more a suitable Dulac function in a part of the interval and the Poincare-Bendixon 
theorem, together with the hyperbolicity of the limit cycle, whenever it exists, for the 
remaining values of m. 

Section [7] deals with the existence of polycycles for the system. Finally in Section [8] 
we glue all the above results to prove Theorem ll.il 
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2. Stability of the origin and proof of Theorem 11.31 

Notice that the origin of (pQ) and ([2]) are nilpotent critical points and there are several 
tools for studying its local stability, see for instance [21 [15l [E|- We will follow the 
approach of [2 [15], based on the polar coordinates introduced by Lyapunov in [T7], to 
study of the stability of degenerate critical points. 

Let u(tp) = Cs(ip) and v(ip) = Sn(y>) be the solutions of the Cauchy problem: 



u = — v 



2 P-\ v = u 2q -\ u(0)= 2 ^/T/p and u(0) = 0, 



where the prime denotes the derivative with respect to ip. 

The Lyapunov generalized polar coordinates are x = r p Cs(ip) and y = r q Sn((p). 
They parameterize the algebraic curves px 2q + qy 2p = r 2pq , that correspond to the level 
sets of above (p,q) — quasi-homogeneous Hamiltonian system. In particular, pCs q {tp) + 
qSn p (tp) = 1, and both functions are smooth Tp^-periodic functions, where 

- r M = 2 P - i /y i /» r _^ r ^ 

— + — 

2p ^ 2q 

and r denotes the Gamma function. The general expression of a differential system in 
these coordinates is: 

x 2q ~ l x + y 2p ~ 1 y ■ pxy — qyx 

y.2pq — l rP+q ' 

In the nilpotent monodromic case, the component 9 does not vanish in a punctured 
neighborhood of the critical point. Hence, system ^ can be written in a neighborhood 
of r = as 

A °° 

^=X>(0y, (5) 

i=\ 

where Ri(0), i > 1 are T-periodic functions. The solution of © that for 9 = passes 
for r = p can be written as the power series 

oo 

r(8,p) = p + J2 u i(8)p\ with «i(0) = 0, (6) 

j=2 

and the functions ui can be computed solving recursive linear differential equations 
obtained plugging © in ([5]) . It is well-known that the stability of the origin is given by 
the first non- vanishing generalized Lyapunov constant Vk := Uk{T). 

To effectively compute some integrals of the above generalized trigonometric functions 
we will use the following result, see [T5] . 



Lemma 2.1. Let Sn and Cs be the (l,q) -trigonometrical functions and let T be their 
period. Then, for i,j £ N, 

(i) J" Sn l (6) Cs J (9) d9 = when either i or j are odd. 

2r(^±i)rfe-) 

(ii) L Sn l (9) Cs J (9) d9 = — -— -, r- when i and j are both even. 

(i„) For q = 2, /„» a» M # = «N»)tf[*)i^oai«lt» 
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Proof of Theorem M.SX By using the transformation (x,y) — > (y,x), system ([2]) becomes 



x = —y + mx 2s+1 , 
y = x 3 — y 2k+1 . 



(7) 



(S) 



We use (HJ), with p = 1 and q = 2, to transform it into 

r = mCs 2s+ \6)r 2s+1 - Sn 2fc+2 (6)r ik+1 , 

9 = r- Cs(0) Sn 2fc+1 (#)r 4fc - 2m Cs 2s+1 (0) Sn(0)r 2s , 

or equivalently, 

dr _ mCs 2s+ \9)r 2s - Sn 2k+2 (9)r Ak 

d9 ~ 1 - Cs(0) Sn 2 ^ 1 ^)^" 1 - 2m Cs 2s+1 (^) Sn^Jr 2 "- 1 ' 

The Taylor series of the right hand-side of ([8]) at the origin has three possibilities: 
(i) When s < 2k, then (|8|) becomes 

^- = mC S 2s+ \8)r 2s + 0(r 4k ). 
do 

Therefore, using the method explained above and Lemma 12.11 we get that its first 

Lyapunov constant is 

V 2s = m Cs 2 ^(6)d9= mV f l s } 4 I . (9) 

Jo 1 (-Z-) 

Then m = is the bifurcation value, and the origin of ^) changes its stability from 
attractor to repeller as m goes from negative values to positive values, 
(ii) Suppose s > 2k, then the Taylor expansion of ([5]) at r = is 

± = -Sn 2k+2 (9)r 4k + 0(r 2s ). 
da 

By using the same method, we obtain that the first Lyapunov constant is 



r(|)rm 



,2fe+2/ i _ 

, 2fc+l . 



V ik = / - Sn 2k + 2 (9) d6 = - 2 \i[ 2 J < 0, (10) 



and the stability of the origin of (|2j) is independent of m and it is an attractor for all m. 
(iii) Finally, when s = 2k we have 

dr 

dd 
Hence the first non-vanishing generalized Lyapunov constant is given by 



% = (m Cs 4fc+ » - Sn 2fc+2 (0)) r 4fc + 0{r* k ~ l ). (11) 

do V / 



V Ak = f (mCs 4fc+4 (0) - Sn 2fc+2 (0)) dd. 



'o 
By using (|9|) with s = 2k and (fT0|) . after some simplifications, we obtain that 

_ 2vr 3 / 2 (m(4A; + 1)!!!! - (2k + 1)!!) 

(r (|)) 2 (4^ + 3)!!!! 

Therefore the origin of ([2]) is attractor for m < (2k + l)!!/(4/c + 1)!!!! and repeller for 
m > (2k + l)!!/(4/c + 1)!!!!, as we wanted to prove. 

In the particular case s = 2k and k = 1, that corresponds to system ([T|), and when 
m = 3/5 we have that V4 = 0. To proof the theorem we continue computing the next 
non-zero Lyapunov constant. For s = 2 and k = 1, equation ([8]) writes as 

± = R 4 (ey + R 7 (9)r 7 + R 10 (6)r 10 + 0(r 13 ), 
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with R A (8) = mCs 8 {8) - Sn 4 (0), 

R 7 {9) = 2m 2 Cs 13 (0) Sn(0) + mCs 9 (6) Sn 3 (0) - 2mCs 5 (8) Sn 5 (0) - Cs(0) Sn 7 (0) 

and 

R w (8) =4m 3 Cs 18 (9) Sn 2 (0) + Am 2 Cs 14 (0) Sn 4 (0) + m(l - 4m) Cs lo (0) Sn 6 (0) 
- 4mCs 6 (0) Sn 8 (0) - Cs 2 (0) Sn lo (0), 

with m = 3/5. Following the procedure explained in this section we obtain that U2 = 
u 3 = 0, 

114(0) = / R A (tp) dip, u 5 = u G = 0, 

Jo 
f e 
u 7 (9) = / (R 7 (i>) + 4R^)u 4 (ijj)) di/>, u 8 = u 9 = 0, 

Jo 



u 10 (0) = / ^io(V') + 7R 7 (i/>)u 4 (ti>) + 4# 4 (VW(V0 + ei^KM) <^- 

Using Lemma ET] and some easy computations we get that V\ = ■ ■ ■ = Vg = 0. Finally, 
it suffices to compute 

Vio = / (i2io(<?) + 7R 7 (8)u A (8) + 4^(0)1x7(0)) d0, 

du\(0) 
because — j- — = 3R A (0)u\(0). Using integration by parts and the expression of u 7 we 

dv 
arrive to 

V ™ = f ( R w(0) + 3«4 (0)u' 7 (0J} dO = j [R lG (9) + 2,u A (8)R 7 (9))d6. (12) 

Notice that applying (iii) and (iv) of Lemma 12.11 we know that 

3 ^8/n „ 4/A \ ,,_6 Sn(0) Cs 5 (0) + 15 Sn(0) Cs(0) + 5 Sn 3 (0) Cs(0) 



U4(0) = J (Jcs 8 (V0-Sn 4 (V>))^ 



35 



Plugging this expression in (|12|) , using several times (i) and (ii) of Lemma 12.11 and the 
properties of the V function we arrive to 

1625 y/ir 

Hence the origin is unstable for m = 3/5. As a consequence, we know that at m = 3/5 
the system has a Hopf-like bifurcation. Therefore the system has at least one limit cycle 
near the origin for m < 3/5. □ 



3. More preliminary results 

This section is a miscellaneous one and it is divided into several short subsections 
containing either some tools that we will use to prove Theorems 11.11 and 11.21 or some 
preliminary results. 
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3.1. Differences between families that are SCFRVF and families that are not. 

On the one hand, if a one-parameter family of differential systems is a SCFRVF, then 
there are many results that allow to control the possible bifurcations; see [9j[22j[23]. One 
of the most useful ones is the so called non-intersection property. It asserts that if 71 and 
72 are limit cycles corresponding to systems with different values of m, then 71 072 = 0. 
Informally, we like to call this property Atila's propertyj because it implies that, if for 
some value of ma limit cycle passes trough a region of the phase plane, this region 
turns out to be forbidden for the periodic orbits that the system could have for any 
other value of the parameter. As a consequence, in this case, the study of 1-parameter 
bifurcation diagrams is much simpler. 

For instance, consider a 1-parameter SCFRVF satisfying the following property: 
(P) For each m £ (mo, mi), it has at most one limit cycle, that we denote by j m . Here, 
if for some m the corresponding system has no limit cycles then 7 m = 0. Moreover, 
assume that U me ( m0jmi )7 m covers a region of the plane where all the periodic orbits of 
the system have to cut. 

Therefore, it holds that: for m E K\ (mo, mi) the system has no periodic orbits. 

The above property has very important practical consequences if we want to deter- 
mine the values mo and mi, that constitute, in many cases, the most difficult ones to 
be obtained to complete the bifurcation diagram. Usually, one of the values, say mo 
corresponds to a Hopf-like bifurcation, and some local analysis allows to obtain it. Then, 
for instance, if for some value of m, say m > mo, the system has no limit cycles then 
m\ < m. The same idea can also be applied to obtain lower bounds of mi. These facts 
simplify a lot the obtention of analytic bounds for the value mi because it suffices to 
deal with concrete systems, with fixed values of m. This approach has been applied 
with success in many works; see for instance [HI Q31 [211 I23JI25]. 

On the other hand, if for a general family of vector fields, we have that the same 
property (P) given above holds, we can say nothing of what happens for m 6 K \ 
(mo, mi). For this reason, when we study system (pQ), we can not ensure the existence of a 
unique value of m for which phase portrait (b) of Figure [1] appears; see also Example l7.1[ 
We remark that system (pQ) is not a SCFRVF with respect to m, and moreover we have 
not been able to transform it into an equivalent one that were a SCFRVF. 

From our point of view, to introduce tools for studying 1-parameter families that are 
not SCFRVF is a challenge for the differential equations community. 

3.2. Global phase portrait. We will draw the phase portraits of system ([1]) on the 
Poincare disc, [3j [23]. Recall that, from the works of Markus [18] and Newmann [20j . 
for knowing a phase portrait it suffices to determine the type of critical points (finite 
and at infinity), the configuration of their separatrices, and the number and character 
of their periodic orbits. 

We start making a study of the critical points at infinity of the Poincare compactifi- 
cation of the system. That is, we will use the transformations (x,y) = (l/z,u/z) and 
(x, y) = (v/z, 1/z), with a suitable change of time to transform system ([1]) into two new 
polynomial systems, one in the (u, z)-plane and another one in the (v, z)-plane respec- 
tively; see [3] for the details. Then, for understanding the behavior of the solutions of (P) 
near infinity it suffices to study the type of critical points of the transformed systems 
which are localized on the line z = 0. These points are precisely the so called critical 
points at infinity of system ([1]). 



Recall that about Atila, King of the Huns, it was said that "the grass never grew on the spot where 
his horse had trod" . 
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Lemma 3.1. By using the transformation (x,y) = (v/z,l/z) and the change of time 
dt/dr = 1/z 4 system (JTJ) is transformed into the system 

z' = —mz + vz 5 , 

where the prime denotes the derivative with respect to r . The origin is the unique critical 
point of (|13p on z = and it is an attracting node. 

The proof of the above result is straightforward. 

Lemma 3.2. By using the transformation (x,y) = (l/z,u/z) and the change of time 
dt/dr = 1/z 4 system ([I]) is transformed into the system 

(u — z 2 )z 2 + u 4 (mu — z 2 ) 
z' = (1 — u 3 )z 

where the prime denotes the derivative with respect to r. The origin is the unique critical 
point of (|14p on z = and it is a repeller. 

Proof. From the expression of (|14p it is clear that the origin is its unique critical point 
on z = 0. For determining its nature we will use the directional blow-up since the linear 
part of the system at this point vanishes identically; see again [3]. 

We apply the z-directional blow-up given by the transformation r = u/z, z = z. 
Performing it, together with the change of time dt/dr = z 3 , system (|14[) is transformed 
into 

r = —1 + mzr 5 , 



'-n „3u3 (14) 



* = l-*V. (15) 

System (|15p has no critical points on z = 0. Then by using the transformation (u, z) = 
(rz, z) we can obtain the phase portrait of system (|15p . Recall that the mapping swaps 
the third and fourth quadrants in the z-directional blow-up. In addition, taking into 
account the change of time dt/dr = z 3 , it follows that the vector field in the third and 
fourth quadrant of the plane (u, z) has the opposite direction to the one obtained in the 
(r, z)-plane. 

Next, we need to perform the it-directional blow-up for knowing the phase portrait 
in such direction. After that, by joining the information about the blow-ups in both 
directions, we will have the phase portrait of system (1141) . 

The it-directional blow-up is given by the transformation u = u, q = z/u, and with 
the change of time dt/dr = u 3 , system (|14p is transformed into 

f(uq 2 - 1) - u 2 (uq 2 - m), .^ 

— muq. 

On q = 0, the origin is the unique critical point of the system, and since the linear part 
of the system at this point vanishes identically we have to use again some directional 
blow-ups. 

Since the lower degree term of qu — iiq is —q(2mu 2 + q 2 ), and it only vanishes on the 
direction q = 0, to study the origin of system (|16p it suffices to consider the n-directional 
blow-up. It is given by the transformation u = u, s = q/u. Doing the change of time 
dt/dr = u, system (fj"6|) becomes 

u = —us 2 (u 3 s 2 — 1) — u(u 3 s 2 — m), , „s 

s = s 3 (u 3 — 1) + 2s(n 3 s 4 — m). 

For s = 0, system ()17p has a unique critical point at the origin. The linearization matrix 
at the origin has eigenvalues m and — 2m. Thus the origin of system (fl"T|) is a saddle. 

Then by using the transformation (u, q) = (u, su) we can obtain the phase portrait of 
system ()16p . Recall that the mapping swaps the second and the third quadrants in the 
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u-directional blow-up. In addition, taking into account the change of time dt/dr = u 
it follows that the vector field in the second and third quadrants of the plane (u, q) has 
the opposite direction to the one in the (u, s)-plane. Once we have the phase portrait 
in the (u, g)-plane, we apply the transformation (u,z) = (u,qu). 

By considering the properties of the blow-up technic and from the analysis of all the 
intermediate phase portraits we obtain that the origin of system (|14[) is a repeller. □ 

Recall that the finite critical points are two hyperbolic saddles at (±m — 1 ' , ±m -1 ' ) 
and a monodromic nilpotent singularity (0,0), whose stability is given in Theorem 11.31 
Finally notice that the vector field is symmetric with respect to the origin. By joining 
to these properties all the information concerning the infinite critical points and using 
the existence and uniqueness results on the number of limit cycles and polycycles given 
in Theorem II. 1[ we obtain the global phase portraits of system ([JJ given in Figure [I] 

3.3. Some Bendixson-Dulac type criteria. Next statement is a Bendixson-Dulac 
type result, that mixes the Bendixson-Dulac Test given in the classical book [31 Thin. 
31] and the one given in |121 Prop. 2.2]. It is adapted to our interests. Similar results 
appear in 131 US |2g] . 

Proposition 3.3 (Bendixson-Dulac Criterion). Let X = (P,Q) be the vector field as- 
sociated to the C 1 -differential system 



x = P(x,y), 
y = Q(x,y), 



;i8) 



and letlA C M 2 be an open region with boundary formed by finitely many algebraic curves. 
Assume that there exists a rational function V(x, y) and k € IR + such that 

M = M {v>k} (x,y) = (VV,X) - kVdiv(X) (19) 

does not change sign in Li and M only vanishes on points, or curves that are not invariant 
by the flow of X. Then, 

(I) If all the connected components oflA\{V = 0} are simple connected then the 

system has neither periodic orbits nor polycycles. 
(II) If all the connected components ofU\{V = 0} are simple connected, except one, 
say U, that is 1-connected, then, either the system has neither periodic orbits 
nor polycycles or it has at most one of them in U. Moreover, when it has a limit 
cycle, it is hyperbolic, is contained in U, and its stability is given by the sign of 
-VM on U. 

Proof. Consider the Dulac function g(x,y) = \V{x,y)\~ l ' k . Then 

d " ( ^ = l p+ | (? + 9( f + f) = <V S ,X) +9 d„(X) 

= -isign(y)|I/|-^ i ((W,X) -jfcVdivm) 
k 

= -^signOOlFr^M^jfc} = -^siga(V)\V\-^M. 

By the hypotheses, M|{y =0 } = (W, X)|{y =0 } does not change sign in U and there is no 
solution contained in {M = 0}. Therefore, neither the periodic orbits nor the polycycles 
of the vector field in U. can intersect {V = 0}. 

For proving (I) we follow the proof of the Bendixson-Dulac Criterion given in [31 
Thm. 31]. Assume, to arrive a contradiction, that the system has a simple closed curve 
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r which is union of trajectories of the vector field. Let C <ZlA the bounded region with 
boundary V. Then, by Stokes Theorem, we have that 



II dlv(gX) = J (gX,n) 



where T is oriented in the suitable way. Note that the right hand-side term in this 
equality is zero because gX is tangent to the curve T and the left one is non-zero by our 
hypothesis. This fact leads to the desired contradiction. 

In case (II), applying a similar argument to the region bounded by two possible simple 
closed curves formed by trajectories of the vector field, we arrive again to a contradiction. 

To end the proof, let us show the hyperbolicity of the possible limit cycle T. Write 
r = J7(t) := (x(t),y(t)),t G [0,T]} C U, where T is its period, and its characteristic 
exponent as h(T) = j Q div (X(7(i))) dt. We need to prove that h(T) ^ and that its 
sign coincides with the sign of —VM on U. We know that 

M (VV,X) 



V V 



kdiv(X). 



Remember that T n {V = 0} = 0. Evaluating this last equality on 7 and integrating 
between and T we obtain that 

T ^(i(t))dt= [ T( ?v^l m)dt _ k f T d[y{X )(^t))dt 

v Jo v Jo 

= ln\V(j(t))\ t ~_ -kh(T) = -kh(T). (20) 

Therefore, the result follows. □ 

Next result is an straightforward consequence of the above proposition. It notices 
that when we construct a suitable Dulac function, the same method provides an effective 
estimation of the basin of attraction of the attracting critical points. 

Corollary 3.4. Assume that we are under the hypotheses of the above theorem and 
moreover that {V(x, y) = 0} has an oval such that it and the bounded region surrounded 
by it, say W, are contained in U. Then, if the differential system has only a critical 
point p in W, and it is an attractor, then W is contained in the basin of attraction of p. 

Observe that when we are under the hypotheses of the above corollary, but we already 
know that the system has a limit cycle in U, and U is simply connected, then, unless the 
set {V(x, y) = 0} reduces to a single point, there is no need to assume that {V(x, y) = 0} 
has an oval. The existence of the oval is already guaranteed by the method itself. 

Sometimes the hypothesis that M does not change sign can be replaced for another 
one, following next remark. 

Remark 3.5. Assume that in Proposition l^.^l we can not ensure that the function M, 
given in (|19|) . keeps sign on the whole domain U. Then, this hypothesis can be changed 
by another one. Define {M = 0}* to be the subset of {M = 0} formed by curves that 
separate the regions {M > 0} and {M < 0}. Then, the new hypothesis is that the 
set {M = 0}* is without contact by the flow of X . Then, in the conclusions of the 
proposition, the connected components ofU\{V = 0} must be replaced by the connected 
components ofU\ ({V = 0} U {M = 0}*) and the same type of conclusions hold. We 
will use this idea in the proof of Proposition \6.1[ 
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3.4. Zeros of 1-parameter families of polynomials. As usual, for a polynomial 

P(x) = a n x n + • • • + a\x + clq, we write A X (P) to denote its discriminant, that is, 

A X (P) = (-1)^ — Res(P(aO,P'(aO), 

where Res(P, P', x) is the resultant of P and P' with respect to x; see [7]. 

By using the same techniques that in |lll Lem. 8.1], it is not difficult to prove the 
following result that will be used in several parts of the paper. 

Lemma 3.6. Let G m (x) = g n {rn)x n + g n -i(m)x n ~ 1 + • • • + gi(m)x + go(m) be a fam- 
ily of real polynomials depending continuously on a real parameter m and set A m = 
(c(m),d(m)) for some continuous functions c(m) and d(m). Suppose that there exists 
an interval Icl such that: 

(i) For some rriQ E I, G mo has exactly r zeros in A mo and all them are simple. 
(ii) For all m £ I, G m (c(m)) ■ G m (d(m)) / 0. 
(iii) For all me I, A x (G rn ) / 0. 

Then for all m G /, G m (x) has also exactly r zeros in A m and all them are simple. 

The idea of the proof consists in looking at the roots of G as continuous functions of 
m. The hypothesis (ii) prevents that some real roots of G m passes, varying m, trough 
the boundary of A m . The hypothesis (iii) forbids, that varying m, appears some multiple 
root of G m . 

Notice that the above result transforms the control of the zeros of a function depending 
on two variables, x and m, into three problems of only one variable, the one of item (i) 
with the variable x and the two remainder ones with the variable m. If the dependence 
on m is also polynomial, and the polynomial has rational coefficients, then these three 
simpler questions can be solved by applying the well-known Sturm method. As we 
will see in the proof of Proposition 15. 2\ this approach can also extended when the one 
variable polynomial has some irrational coefficients. 

3.5. Transformation into an Abel equation. System (pQ) can be seen as the sum 
of two quasi- homogeneous vector fields, see [6]. It is known that in many cases these 
systems can be transformed into Abel equations. We get: 

Proposition 3.7. The periodic orbits of system (HJ correspond to positive T-periodic 
solutions of the Abel equation 

± = a{ e )p 3 + mp 2 , (21) 

where 

a{6) = 3 Cs(0) Sn(0) (2m Cs 4 (0) + Sn 2 (#)) (m Cs 8 (0) - Sn 4 (#)) 
and 

P{6) = 5m Cs 8 {9) - 4Sn 4 (0) + (3 - 10m) Cs 4 (#) Sn 2 (0), 
being Sn and Cs the functions introduced in Subsection^ and T their period. 

Proof. The result follows by applying the Cherkas transformation 

P ~ 1 - r 3 Sn(0) Cs(0) (Sn 2 (#) + 2m Cs 4 (0)) ' 

to the expression of system (pQ) in the quasi-homogeneous polar coordinates introduced 
in Section [2j It is used that the periodic orbits of the system do not intersect the curve 
= 0, and therefore the above transformation is well-defined over them, see [6]. □ 
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Using the above expression it is not difficult to reproduce the proof of the existence of 
the Hopf-like bifurcation given in Subsection [2j Unfortunately, although expression ()2ip 
looks quite simple, the results about the number of limit cycles of Abel equations that 
we know are not applicable to ()2ip . 



4. NON-EXISTENCE OF LIMIT CYCLES FOR m G (0, 9/25) U (3/5, Oo) 

In this section we prove the non-existence results of periodic orbits already given in 
[10] and extend them to the non-existence of polycycles. Our proof is different and based 
on the Bendixson-Dulac theorem and other classical tools. We study separately each 
interval. 

Proposition 4.1. Form G (0,9/25], system ([1]) has neither periodic orbits nor polycy- 
cles. 

Proof. Recall that for m G (0, 9/25] the origin is attractor. Therefore if we prove that 
any periodic orbit V of the system is also attractor we will have proved that the system 
has not periodic orbits. In order to prove the stability of the limit cycle we need to 
compute J" div(X(7(t))) dt, where ^{t) := (x(t),y(t)) is the time parametrization of 
r, and T = T(T) its period. 

From equation (|19p . for any function V such that {V(x, y) = 0} n T = 0, we have 



-kV 
Hence, 



div (x( 7 (t))) dt=~ rs/ ,:: dt, 

o Jo kV{j{t)) 

where we have followed similar computations that in (|2U|) . Then the stability of T is 
given by the sign of —MV. If we show that for m G (0, 9/25] there exist a non-negative 
V and k G M + , such that its corresponding M is non-negative, then we will have proven 
that the limit cycle is hyperbolic and attractor. 

By considering V(x, y) = 2x 2 + y 4 and k = 2/3 equation (|19|) becomes 

M {v,§} = I (( 3 - 10m ) x2 + m v 4 ) v*> 

which clearly is non-negative on M? for m G (0,3/10]. 

If we use the same V(x, y) that in previous case, but k = K(m) = 8(llm + it!)/(10m-|- 
3) 2 , with R = ^/m{l — 4m) (25m — 9) then we have 

(m + i?)(llm + R)\ 1/2 2 2(3-10m) fm(llm + R)\ 1/2 4 \ 



M{V ^ (m)} - [JTIQ^ { m ) X + 3 + 10m v (m + R) ) V , 

hence M^y K ^ m ^ is non-negative on M 2 for m G (1/4,9/25]. 

Therefore system (JTJ has no limit cycles for m G (0, 9/25] as we wanted to show. 

To prove the non-existence of polycycles for m G (0, 9/25) we use a different approach. 
Following [24] , we can associate to each polycycle T, with k hyperbolic saddles at its 
corners, the number p(T) = Y\ i=1 bi/ai, where — ai < < 6j, i = 1, . . . k, are the eigen- 
values at the saddles. Then, T is stable (respectively, unstable) if p(T) < 1 (respectively, 
p(r) > 1). In our case 

( 5<Jm - 3 + ^25771 + 18^ + 9) 

" (r) = ^ 48^ "■ 

Then, easy computations show that the polycycle is an attractor if m < 9/25 and 
a repeller if m > 9/25. Assume, to arrive to a contradiction, that for m < 9/25 the 
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polycycle exists. Then both, the polycycle and the origin, would be attractors. Applying 
the Poincare-Bendixson Theorem we could ensure that the system would have at least 
one periodic orbit between them. This result is in contradiction with the first part of 
the proof, where the non-existence of periodic orbits is established. 

It only remains to show that for m = 9/25 the polycycle neither exists. To prove this 
fact we could study the stability of the polycycle showing that if it exists it would be 
attractor, arriving again to a contradiction. Nevertheless it is easier to apply Proposi- 
tion [33] with the V and k = K{9/25) used to prove the non-existence of periodic orbits. 
Indeed, this later approach, taking the corresponding V and k, could also be used for 
all values of m £ (0, 9/25], but we have preferred to include a proof based on the study 
of the stability of the limit cycle and the polycycle. □ 

Lemma 4.2. Let X be the vector field associated to system ([1]). 

(i) If we take k = 1/3 and Vi(x,y) = g (y) + g x (y)x where g (y) = g[(y) and g x (y) 
a solution of the second order linear ordinary differential equation 

- g'l(y) + my 5 g[(y) - ^my A gi (y) = 0, (22) 

then (|19[) reduces to the function 

Mi := M {Vu i } (x,y) = ^y 3 (3my 2 g'((y) - 5myg' 1 (y) + 3gi(y)) . (23) 

(ii) If we take k = 2/3 and V 2 {x,y) = go{y) + 9i(y)x + <?2 {y)x 2 , with 

9i(y) = 92(y), 

g (y) = (l/2)g' 2 \y) - (l/2)my 5 g' 2 (y) + (5/3)m/ 52 (y), (24) 

then (|19p becomes 

-L /// , \ 3 a ii , \ 5 a i i \ 2, , 



M 2 := M {V2 2 } (x,y) = ( - -g 2 "{y) + -my 5 g' 2 \y) - -my 4 g 2 {y) + -(3 - 10m)y 3 52 (y))x 

))• 
(25) 



+ -3-y 3 (9my 2 g' 2 "(y) - m(30 + 9my 6 )yg 2 \y) + 3(6 + 5m 2 y 6 )<7 2 (y) + 20m 2 y 5 g 2 {y)) . 
lo 



Proof (i) If Vi(x, y) = g (y) + g 1 (y) x and k = 1/3, then 
M 1 =(VV 1 ,X)-^dw(X)V 1 

= {go(y) - 9i{y))x 2 + y- g' {y) + my 5 'g[(y) - -my 4 gi(y)jx 

+ ^V 3 {^mg' (y)y 2 - bmg (y)y + 3gi(y)) . 

By choosing go(y) = g'\{y) the coefficient of x 2 in M\ vanishes, and we obtain 

Mi = f -g'l (y) + my 5 g[(y) - -™y 4 .gi(y) J x + -i 

Finally, if g\{y) is a solution of (|22p we get (|23p . 
(ii) If k = 2/3 and V 2 {x,y) = g {y) + g\{y)x + g 2 {y)x 2 , then 

M 2 =(V^ 2 ,X)-^div(X)T/ 2 

=(9i(y) - g' 2 {y))x 3 + (rny 5 g' 2 {y) - —my 4 g 2 (y) - g[(y) + 2 5o (y))x 2 + (2y 3 g 2 {y) 
+ my 5 g[(y) - —my^giiy) - g' (y)jx + -y 3 (3gi(y) + 3my 2 g' (y) - Wmyg (y)). 



Ah = f -g"{y) + m y 5 9i(y) - ^ m y 4 9i(y) J x + -y 3 (3mg"(y)y 2 - 5myg[(y) + 3gi{y)) . 
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By choosing gi(y) = g' 2 (y) and g (y) = (l/2)g%(y) - (l/2)my 5 g' 2 (y) + (5/3)my 4 g 2 (y) the 
coefficients of x 2 and x 3 in M 2 vanish. Then we have (|25p. D 

Remark 4.3. Notice that if 52(2/) is a solution of the linear ordinary differential equation 
~ \&'{v) + l^y 5 92(y) - lrny A g 2 (y) + |(3 - I0m)y 3 g 2 (y) = 0, (26) 



then (|19p reduces to a function depending only of the variable y. 

Proposition 4.4. Form E [3/5, 00), system ([1]) has neither periodic orbits nor poly cy- 
cles. 

Proof. We want to apply Proposition 13.31 taking k = 1/3 and V\(x,y) = go(y) + gi(y)x, 
with go and g\ as in (i) of Lemma 14.21 Applying the transformation z = my e /6, 
equation (f22|) becomes 



**?(*)+ (§ -*)tfl(*) + ^ 



<&(*) = 0, 

which is a Kummer equation, see [U pp. 504]. A particular solution of this equation is 

j=o ^eh J- 
where (a),- := a(a + l)(a + 2) ■ ■ ■ (a + j — 1) and (a)o = 1. Therefore we consider 



OG 



mx 1 / 6 r— v (— g)j fm\j y 



9l(y) = ( 1 r) vYl 



which is convergent on the whole K and satisfies (|22|) . Its derivatives are 

m\i/6^(-i) j /my' y 6j ' 



j=0 ^6 A? J 



my/C (-g)j (m\J ^y 6j 1 



7=0 V 6A? J 



Replacing the above functions in (|23p we obtain 
3 — 5m \ /m\ V 6 4 



/ <J — am \ fm\ 



Since (— n)j is negative for all j, it follows that M\ < for m > 3/5, and vanishes only 
on y = 0. Therefore the result follows by applying Proposition 13.31 □ 

5. Uniqueness and hyperbolicity of the limit cycle for m G (1/2,3/5) 

In this section we prove that for m 6 (1/2, 3/5), system ([TJ has at most one limit cycle 
or one polycycle and both never coexist. Moreover, we show that when the limit cycle 
exists, it is hyperbolic. The uniqueness of the limit cycle was already proved in [101 . Our 
approach is different and, like in the previous section, it is based on the construction of 
a suitable Dulac function. This section ends with the proof of Proposition 11.21 
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Lemma 5.1. Let S be the open set bounded by the lines x = ±m -1 / 4 and y = ±m -1 ' 4 
and let O be the connected component containing the origin and bounded by the above 
four straight lines and the hyperbola xy + 1 = 0, see Figure^ The following holds: 

(i) The vector field X associated to system (JTJ) is transversal to the boundary dS of 

the square S except at the two saddle critical points of system (JJ) . 
(ii) If system (pQ) has a periodic orbit or a poly cycle, it must be contained in Q C S. 





(a) (b) 

Figure 3. Regions O and S. 



Proof, (i). Consider the function f(x,y) 



m 



' . It is not difficult to see that 



(Vf,X) restricted to x 



m 



-1/4 



has the expression y 



m 



' which is negative 



for y £ (— m -1 ' 4 , m -1 ' 4 ). Analogously, we can see that the direction of X along dS is 
as showed in Figure [3] (a) . 

(ii) . It is well-known that the sum of the indices of all the singularities surrounded by 
a periodic orbit, or a polycycle is one. Recall that the indices of the saddle points are —1 
and the index of a monodromic point is +1. Hence, if a periodic orbit or a polycycle T 
exist they must surround only the origin. Moreover, by statement (i), V cannot intersect 
dS. Finally, a simple computation shows that (V(xy + l),X) restricted to xy + 1 = is 
(1 — m)/x 4 , which implies that X is transversal to xy + 1 = 0. Hence X is transversal 
to 90 and the lemma follows. □ 

Proposition 5.2. For m £ [1/2,3/5), system ([TJ has at most one limit cycle and one 
polycycle and both never coexist. Moreover, when the limit cycle exists it is hyperbolic 
and repeller. 



Proof. Following statement (ii) of Lemma l4.2l we take k = 2/3 and a function Vi(x, y) = 
9o(y) + 9i(y) x + 92{y)% 2 adequate to apply Proposition 13.31 for proving the uniqueness 
of the limit cycles or polycycles for system (JJJ. 

We will take g^iy) as a truncated Taylor series at the origin of a suitable solution of 
(|26|) such that the curve {V2 = 0} has an oval surrounding the origin, and that M^ does 
not change sign in O. These two properties will imply the result. 

The general solution of ([26]) is the linear combination of generalized hyper geometric 
functions 

1/6 j 



9z(y) = c o^2 



"M)j(^ ( m ))j /my y 



3=0 



+ C 2 y 2 J2 



{3)3 (g)j 

(V> + (m)) • (V>~ 



{m)) j 



6/7 V3J- 



,&)j 



+Ciy S (i), (i), wit 



ly 



(\i 



(27) 
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where ^(m) = ±A(m) - 2/9, ip ± {m) = ±A(m) - 1/18, ^(m) = ±A(m) + 1/9, with 
A(m) = y/(Um - 3)/m/9. 

We look for an even solution, so we take C\ = 0. As we will consider Cq ^ 0, it is 
not restrictive to choose Cq = 1. Finally, the constant C*2 = —(3/5 — m) 2 ' 3 is fixed 
after some previous simulations and taking into account that we already know that at 
m = 3/5 there is a Hopf-like bifurcation. 

Once we have fixed the above constants, we calculate the Taylor polynomial of de- 
gree 12 of (/2 at y = 0, 712(52), obtaining 

T12M?/)) =7^(3 - 10m)(3 + 35m)y 12 - JL(75 - 125m) 2 / 3 (3 - 13m)y 8 

+ -U3 - 10m)y 6 - 4(75 - I25m) 2 / 3 y 2 + 1. (28) 

9U 25 

So, in (ii) Lemma 14.21 we fix 52 as 712(52(2/))- Then the corresponding 50 and 51 are 
given by (|24p . Thus, M2 is of the form M2 = ((j)(y)x + ■0(y))y 4 where 

0(y) =777^ (7^ (3 - 10 m) (242 m + 3) (35 m + 3) y 11 

+ (75 - 125 m) 2/3 (86 m + 3) (13 m - 3) y 7 

947 1 S 

W,( y ) = m 2 (3 - 10 m) (35 m + 3) y 16 m 2 (75 - 125 m) 2/3 (13 m - 3) y 12 

^ vy; 400950 v M ;y 4050 v ; y Jy 

+ (3 - 10 m) (550 m 2 + 145 m + 3) y 10 

7425 v y 

+ (75 - 125m) 2/3 (196 m 2 - 45 m - 9) y 6 

4725 v ' y ' 

1 9 

+ — (3 - 5 m) y 4 (75 - 125 m) 2/3 (3 - 5 m) . 

15 75 

The proposition follows if we prove that M2 does not change sign on the region Vt. In 
fact, it is sufficient to prove that M := A^/y 4 does not change sign on 0. 

The idea is to show that {M = 0} does not intersect Vt. Since M is linear in the 
variable x, {M = 0} cannot have ovals inside £1. If {M = 0} has a component in f2, this 
component would have to cross d£l by continuity of the function. Then, it suffices to see 
that {M = 0} does not intersect d£l. Moreover, as M satisfies M(x,y) = M(—x,—y), 
it is sufficient to study M on half of dil. To deal only with polynomials we introduce 
the new variables n = ^fra and s = (75 — 125m) 2 ' 3 . Notice that s 3 = (75 — 125n 4 ) 2 . 

We split the half of the boundary of £1 in four pieces: 

• The segment 71 = {(x, 1/n) : — n < x < 1/n}, 

• The segment 72 = {(1/n, y) : — n <y < 1/n}, 

• The piece of hyperbola 73 = {(x, —1/x) : n < x < 1/n}, 

• The corners 74 = {(1/n, 1/n), (1/n, — n), (n, —1/n)} 

and we have to prove that {M = 0} n 7$ = for each i = 1, 2, 3, 4. 
These facts can be seen proving that for n G I := [yl/2, y3/5), 

• Qi(x,n, s) := M(x, 1/n) 7^ 0, for x E (— n, 1/n). 

• Q2(y,n,s) := M(l/n,y) / 0, for y G (-n, 1/n). 

• Q'iix, n i -s) := M(x, —1/x) 7^ 0, for x S (n, 1/n). 

• M(l/n, 1/n) • Af (1/n, -n) • M(n, -1/n) / 0. 

Lemma 13.61 with r = 0, is a convenient tool to prove the first three items. The proof 
of the last item is an straightforward consequence of Sturm method. 



18 J. D. GARCIA-SALDANA, A. GASULL, AND H. GIACOMINI 

We will give the details of the proof that Q2(y, n, s) 7^ 0, which is the most elaborate 
case. The remainder two cases follow similarly. 

Writing Q(y,n,s) := 280G650nQ2(y,n,s) we get that 

Q(y, n, s) =1729n 9 (35n 4 + 3)(10n 4 - 3)y 16 - 9009n 9 s(13n 4 - 3)y 12 

- 21(10n 4 - 3)(242n 4 + 3)(35n 4 + 3)y n 

- 378n(10n 4 - 3)(550n 8 + 145n 4 + 3)y 10 + 297s(13n 4 - 3)(86n 4 + 3)y 7 

+ 1188ns(196n 8 - 45ra 4 - 9)y 6 - 187110n(5n 4 - 3)y 4 + 74844ns(5n 4 - 3). 

Looking at Lemma 13.61 with r = 0, it suffices to prove the following three facts: 

(i) When n = ^/lj2 G I, Q(y, n, s) / for y G (-n, 1/n). 
(ii) For n£ I, A y Q(y, n, s) / 0. 
(iii) For n £ I, Q(—n, n, s) ■ Q(l/n, n, s) ^ 0. 

Since the polynomial has no rational coefficients the proof of item (i) requires some 
special tricks. Notice that when n = yl/2 then s = 5\/10/2. Hence, 

R(y) := Q(y, ± |^io) =^vW 6 - H|ii^lo^ - imiuy^ 

Ar - in 239085 3/ — 7 51975 4r - 3/ — fi 
- 80514 v^y 10 + — - — \Zl0y 7 + — — vWlOr 6 



93555 ^ 93555^^; 

2 2 

We will prove that the above polynomial has no real roots in [— 1, 12/10] D (—w, 1/ra). 
The Sturm method gives polynomials with huge coefficients and our computers have 
problems to deal with them. We use a different approach. We know, that 

3002 4/ - 37 28 3/ — 265 

n := < V8 < — =: n, s := — < V 10 < =: s, 

1785 22 ' - 13 123 

where these four rational approximations are obtained computing the continuous fraction 
expansion of both irrational numbers. If we construct the polynomial, with rational 
coefficients, 

R+{y) = 70889_ yl6 _ 315315^ 12 _ m7Uy u _ 805Ujiy io 

239085_ - 51975 - 93555_ 4 93555 

H 2 — sy H 2 — HSy H ~ ny 2~ --' 

it is clear that for y > 0, R(y) < R + (y). In fact, 

, , . 2622893 1fi 2427117 12 „„ , „ 11509668 in 

R+(y) = y 16 y 12 - 106764?/ 11 y 10 

\yj 176 y 68 y y 85 y 

21119175 7 15442875 fi 314685 4 37446948 
82 y 164 y 4 tf 221 

and, now, using the Sturm method it is quite easy to prove that R + (y) < for y G 
[0, 12/10]. Hence, in this interval, R(y) < R + (y) < 0, as we wanted to prove. 
To study the values of y < we construct a similar upper bound, 

BT(v) =™^ny 1B - ^^nsy 12 - 106764y n - 80514ny 10 

239085 7 51975 _ - 93555 _ 4 93555 

H ^ — - y H — ns V H — n 2/ g - --' 

and applying the same method the result follows. 
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To prove (ii) we compute 

A y Q{y, n, s) = n 4 V(5n 4 - 3) 5 (35?i 4 + 3) 3 (10n 4 - 3) 3 P 2 58(n, s), 

where -P25s( n ; s ) is a polynomial in n and s of degree 258. Clearly, the roots of the 
first five factors of the above discriminant are no relevant for our problem because the 
corresponding n is not in /. To study whether P 2 58( n , s ) vanishes or not we compute 

Res(P 25 8(n, s), (75 - 125n 4 ) 2 - s\ s) = (5n 4 - 3) 24 P 3 9o(™ 2 ), 

where i-3go(n 2 ) is a polynomial of degree 390 in n 2 . Applying again the Sturm method 
we get that P^q{ji 2 ) has no significant roots for our study. Finally, the numerator of 
Q(—n, n, s) ■ Q(l/n, n, s) is a polynomial in n and s of degree 49. Using the same trick 
as above we prove item (iii). In this case the polynomial we have to deal with has degree 
152 in n. 

Therefore {M = 0} n dQ = and as a consequence {M = 0} n = 0. 

Finally, it is not difficult to see, because V is quadratic in x, that the set {V(x, y) = 0} 
has exactly one oval surrounding the origin. Hence, the proposition follows. □ 



Proof of Proposition \1.S\ Notice that the function V used in the proof of Proposition 
coincides with the function V(x, y, m) of the statement of the proposition. Taking 
k = 2/3 we are also under the hypotheses of Corollary 13.41 Therefore the set lA m is 
contained in Wq, as we wanted to prove. □ 

We remak that following similar ideas that in the above proof we can construct bigger 
sets contained in Wq. For a given m, let us denote by Ti(g 2 (x; C 2 )) the Taylor polynomial 
of degree £ at x = 0, of the function (|27p with Cq = 1, C\ = 0. Then for each £ G N and 
C2 £ R we can take this function as a new seed 52 for constructing the corresponding 
V as in (ii) of Lemma [4.21 Then checking that the oval contained in {V = 0} is crossed 
inwards by the flow of the system, the result follows for the function V constructed with 
these £ and C2. 

6. NON-EXISTENCE OF LIMIT CYCLES AND POLYCYCLES FOR m G (9/25, 0.547] 

This section contains new non-existence results for system (pQ). We split the interval 
into the subintervals (9/25, 1/2) and [1/2, 0.547]. Recall that our numerical study shows 
that the system has no limit cycles for m < 0.56011 ... As m becomes closer to this 
bifurcation value the proof of non-existence of periodic orbits and polycycles becomes 
harder. 

Proposition 6.1. Form £ (9/25,1/2), system ([1]) has neither limit cycles nor polycy- 
cles. 

Proof. We would like to apply Proposition 13.31 To this end we will follow similar steps 
that in the proof of Proposition 15.21 but with a function V such that the set {V = 0} 
has no oval in Q. Recall that f2 is the domain introduced in Lemma I5.1L where the limit 
cycles and the polycycles must lay. We take V = V2(x,y) = go(y) + 9i(y)x + Q2{y)x 2 
with gi(y) = g' 2 (y), g = (l/2)g%(y) - (l/2)my 5 g' 2 (y) + (5/3)my 4 g 2 (y). Now we consider 
52 (y) = ao + a 2V 2 + «4y 4 + a%y^ + a$y 8 , with coefficients to be determined. From 
statement (ii) of Lemma r4.2l it follows that the corresponding M 2 is a polynomial function 
in x of the form M 2 = <j)(y)x+ip(y) where 4>(y) and ip(y) are polynomials in the variable y 
whose coefficients depend on a 2 j, j = 0, 1, . . . , 4. In order to simplify the computations, 
we change the parameter m by n 4 to transform V into a polynomial in the variables x, 
y, and n. Since m G (9/25, 1/2) we can restrict our study to n G (0.77,0.844). 

We consider the values of 04, a^ and os such that 4>{y) has a zero at y = of multi- 
plicity nine, we choose the value of a 2 by imposing that M 2 vanishes at the two saddle 
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points of the system and, finally, we use the freedom of changing 52(2/) by \g 2 (y), for 
any 7^ A G M, to remove all the denominators. We obtain that 

g 2 (y) =270(9 + 51n 2 + 213n 4 + 535n 6 + 924n 8 + 756n 10 ) 

- 756n 2 (9 + 42n 2 + 105n 4 + 130n 6 )y 2 

+ 3(3 - 10n 4 )(9 + 51n 2 + 213n 4 + 535n 6 + 924n 8 + 756n 10 )y 6 

- 3n 2 (3 - 13n 4 )(9 + 42n 2 + 105n 4 + 130n 6 )y 8 . 

The corresponding M 2 is of the form 

M 2 (x, y) = ^y 4 (</>(y)x + ij>(v)) =: ^y 4 M(x, y), (29) 

where 

(j)( y ) =3(3 - 10n 4 )(3 + 35n 4 )(9 + 51n 2 + 213n 4 + 535n 6 + 924n 8 + 756n 10 )y 5 

- 3n 2 (3 - 13n 4 )(3 + 86n 4 )(9 + 42n 2 + 105n 4 + 130n 6 )y 7 , 
if}( y ) = - 756n 2 (3 - 5n 4 )(9 + 42n 2 + 105n 4 + 130n 6 ) 

+ 27(3 - 5n 4 )(9 + 51n 2 + 213n 4 + 535n 6 + 924n 8 + 756n 10 )y 4 

- 12n 2 (9 + 42n 2 + 105n 4 + 130n 6 )(9 + 45n 4 - 196n 8 )y 6 

- 40n 8 (3 - 10n 4 )(9 + 51n 2 + 213n 4 + 535n 6 + 924n 8 + 756n 10 )y 10 
+ 91n 10 (3 - 13n 4 )(9 + 42n 2 + 105n 4 + 130n 6 )y 12 . 

Recall that the main hypothesis in Proposition 13.31 is that M does not change on il. 
As we will see, this happens only for n G J := (0.77, n] where n ~ 0.8045592 will be 
precisely defined afterwards. When n £ K := (n, 0.844) the result will be a consequence 
of the variation of Proposition 13.31 described in Remark 13.51 

For n G J, following similar steps that in the proof of Proposition 15.21 we divide half 
of the boundary of £1 in five pieces: 

• The segment 71 = {(x, 1/n) : — n < x < 1/n}, 

• The segment 72 = {(l/n,y) : — n < y < 1/n}, 

• The piece of hyperbola 73 = {(x, — 1/x) : n < x < 1/n}, 

• The corners 74 = {(1/n, —n), (n, —1/n)}, 

• The corner 75 = {(1/n, 1/n)} 

and we will prove that {M = 0} n ji = for each i = 1,2,3,4 and that although 
(1/n, 1/n) G d£l, the set {M = 0} does not enter in £1. From these results we will have 
proved that M does not change sign on Q and, as a consequence, the proposition will 
follow for n G J. 

To prove the fifth assertion it suffices to study the function M in a neighborhood of 
the point (1/n, 1/n) G d£l. By the construction of M, it holds that M(l/n, 1/n) = 0. 
By computing the partial derivatives of M at this point we obtain which is the tangent 
vector of the curve at (1/n, 1/n). Then, it is easy to see that when n G J , in a punctured 
neighborhood W of (1/n, 1/n), it holds that W Pi {M = 0} n O = 0. In fact, n G dJ is 
a solution of the equation 

(dM(x,y) \ 

num = 0, 

\ OX (x,y)=(l/n,l/n)/ 
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where num(-) denotes the numerator of the rational function. Moreover, 

M (x, -) = - 9(nX 7 1) f88200n 16 + 107800n 14 - 4930n 12 
V nJ n 4 V 

- 37380n 10 - 15855n 8 - 2736n 6 + 576n 4 + 108n 2 - 27) (30) 



and n is also the positive root of the polynomial in n appearing in the right-hand side 
of the above formula. Notice that when n = n, the straight line {y = 1/n} is a subset of 
{M = 0}. This fact is the reason for which this approach only works for n E J = (0.77, n]. 
Let us prove the remainder four assertions. As in the proof of Proposition 15.21 they 
follow by showing that when n E J, 

• Ri(x,n) := num(M(x, 1/n)) ^ 0, for x E (— n, 1/n). 

• R 2 (y,n) : = num(M(l/n, y)) / 0, for y E (— n, 1/n), 

• Rs(x,n) := num(M(x, — l/x)) ^ 0, for x E (n, 1/n). 

• M(l/n, -n) ■ M(n, -1/n) / 0. 

That Ri has no zeros in J, is an straightforward consequence of (|30p . 

To study R 2 and R3 we will use Lemma 15. 11 We start computing the discriminants, 

S 2 (n) = A y {R 2 (y,n)), 5 3 (n) = A x (R 3 (x,n)), 

and analyze whether they vanish or not on J. Using the Sturm method we get that 
on J, S 2 vanishes only at one value n 2 ~ 0.8040188 and S3 also vanishes only at one 
value 113 w 0.8045576. The root n 2 of S 2 forces to split the study of R 2 (y, n) in the three 
subcases: n E (0.77, n 2 ), n = n 2 and n E (ri2,n\. Doing the same type of computations 
and reasoning as in the previous section we can prove all the above assertion when 
n ^ n 2 . The case n = n 2 follows by continuity arguments, because it can be seen that 
in this situation R 2 has a real multiple root but it is not in (— n 2 , l/n 2 ). The study of 
-R3 is similar to the one of R 2 and we omit it. We also get that R3 neither vanishes on 
(n, 1/n). 

That for n 7^ n, M(l/n,—n) ■ M(n, — 1/n) ^ is once more a consequence of the 
Sturm method. 

Therefore, when n E J, we are under the hypotheses of Proposition 13.31 and we will 
know that the system has no limit cycles once we have proved that the set {V = 0} has 
no ovals. We defer the proof of this fact until we have considered the case n E K = 
(n, 0.844). 

When n E K, we know that {M = O}n(7^0 and we are no more under the 
hypotheses of Proposition 13.31 Let us see that we can apply the ideas of Remark 13.51 
To this end we have to prove that {M 2 = 0}* n O is without contact for the flow of X. 
Note that {M 2 = 0}* = {M = 0}*. 

We need to show that M = (VM, X) does not vanish on {M = 0}* n 9,. We study 
the common points of {M = 0} and {M = 0} and prove that they are not in 9. First, 
we compute 

M(x,y) = (VM{x,y),X{x,y)) =: y 3 N(x,y), 
and we remove the factor y 3 . We do not care about the points on {y = 0} because 

M(x, 0) = 756n 2 (5n 4 - 3) (9 + 42n 2 + 105n 4 + 130n 6 ) / 0, 

for nE (0,0.88]. 

The resultant Res(M, N, x) factorizes as 

Res(M,iV,a;) = y 2 (n 2 y 2 - l)(P„, 2 (y))(P„,3 4 (y)), 
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where P n ,2(y) an d Pn,34(y) are polynomials in the variable y with respective degrees 2 
and 34 and whose coefficients are polynomial functions with rational coefficients in the 
variable n. 

Clearly, (n 2 y 2 — 1) does not vanish on — 1/n < y < 1/n. By using once more 
Lemma 13. 61 it is not difficult to prove that P n ,2{y) does not vanish either on — 1/n < y < 
1/n, for n G (n, 0.844). Hence we will focus on the factor P n ,34(y)- 

We will use again Lemma I3U1 By using the Sturm method we get that A y (P n ^4(y)) 
has no zeros in the interval K. In fact one zero is n G dK and another one is n* ~ 
0.8445 G" J and this is the reason for which we can only prove the result until n = 
0.844 < n*. By using Sturm method, it can be shown that -P n .34( — 1/n) ■ P n ,34:(l/n) 7^ 
for all n G K and, for instance, for n = hq = 83/100 G K, the polynomial -P no ,34(y) has 
exactly two (simple) zeros in —1/no < y < 1/no- Then, Lemma 13.61 with r = 2, implies 
that P n ;a{y) has exactly two (simple) zeros in —1/n < y < 1/n, for all n G K. We call 
them y = yi(n),i = 1,2 and they are continuous function of n. Therefore, we need to 
prove that the corresponding points in {M = 0} n {N = 0} are outside of Vt. 

Notice that because of the expression of M, given in (|29p . the points in {M = 0} are 
on the curve T = { (— -rrr, y) '■ y G M\{0}}. Moreover it can be easily seen that <j>{y) / 
on the region that we are considering. Therefore the points in {M = 0} n {N = 0} are 
given by the two continuous curves 

For a fixed n G K it is not difficult to prove that the points in 7$, i = 1, 2 are outside of 
fi. If for some n G K there was a point inside £1, by continuity it would be at least one 
point in one of the pieces of boundary of fi formed by the straight line {x — 1/n = 0} 
and the hyperbola {xy + 1 = 0}. To prove that such a point does not exist we compute 
the following two resultants 

Res (num ( - — — J, P n ,3i(y),yj = Pi25o(n), 

Res (num(- y—pr + 1) , Pn,u(v) , VJ = Pl26o(ji), 

where Pg(n) are given polynomials with rational coefficients and degree £. Both polyno- 
mials factorize in several factors and, using once more the Sturm method, we can easily 
prove that do not vanish on K. Hence, {M = 0} n {N = 0} (1 O = which implies that 
{M = 0} n £1 is without contact by the flow of X, as we wanted to prove. 

Since M is linear in the variable x, {M = 0} cannot have ovals. Therefore, by 
Remark 13.51 to end the proof we need to show that the set {V = 0} has no ovals 
either in Q. We claim that the set {V = 0} n Q is without contact by the flow of the 
system. If this happens and {V = 0} had an oval then it would be without contact. 
Then by the Poincare-Bendixson Theorem it should surround the origin. However, by 
considering the straight line passing through the origin y = 9x/10 it is easy to prove, by 
using again Lemma 13.61 that the function V(x, 9x/10) does not vanish on the interval 
-1/n < x < 1/n for all n G (0.77, 0.844). Thus, {V = 0} n {y - 9x/10 = 0} = 0. Hence, 
V has no ovals inside O, as we wanted to see and the proposition follows by using all the 
above results and the reasoning explained in Remark 13. 5[ 

To prove the above claim, it suffices to see that {M = 0} n {V = 0} n £1 = 0. This is 
because precisely, ML„„, = V. 

Recall that when n G J = (0.77, n] then {M = 0} n O = and so the result follows. 
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Let us consider the case n G K = (n, 0.844). To study if {V = 0} and {M = 0} 
intersect, we compute the resultant of M and V with respect to x. We have 

Res(V,M,x) = (n 2 y 2 - l)P n>ao (v), 

where P n ,3o(y) is a polynomial of degree 30 and whose coefficients are polynomial func- 
tions in the variable n with rational coefficients. We want to prove that Res(V, M, x) 
does not vanish on the interval —1/n < y < 1/n for n G K. It suffices to study P n ,3o(y)- 
We will use once more Lemma 13.61 

The polynomial P n ,3o(— l/ n ) ■ -Pn,3o(l/ n ) has no real roots when n G K. Moreover 
hypothesis (i) of Lemma 13.61 holds with r = (no real roots) by considering for instance 
no = 82/100. To see that condition (iii) of the lemma holds, we compute A y (P n> so(y)). It 
is a polynomial of degree 2728 in the variable n which factorizes in several factors, being 
the largest one of degree 594. From this decomposition we can prove that A y (P n ^o(y)) 
has no zeros for n G K. Therefore, by Lemma 13.61 we conclude that P n ,3o(y) does not 
vanish on the whole interval —1/n < y < 1/n for n G K, and the claim follows. □ 

Proposition 6.2. For m G [0.5,0.547], system (pQ) has neither limit cycles nor polycy- 
cles. 

Proof. We will construct a positive invariant region 1Z having the two saddle points in its 
boundary. As we will see, the proposition follows once we have constructed this region, 
simply by using the uniqueness and hyperbolicity of the limit cycle, whenever it exits. 
We remark that in this proof we will not use the Bendixson-Dulac theorem. 

Assume that such a positive invariant region 1Z exits. By the Index theory, if the 
system had a limit cycle, it should surround only the origin. By Proposition 15.21 we 
already know that for n G [0.5,0.6) D L := [0.5,0.547], the limit cycle would be unique, 
hyperbolic and repeller. By the Bendixson-Poincare Theorem the above facts force 
the existence of another limit cycle and so a contradiction. It is straightforward that 
the existence of this positive invariant region is not compatible with the existence of a 
polycycle connecting both saddle points. 

To construct 1Z we consider a function V2(x,y) = go(y) + 9i{y)x + g2(y)x 2 , with go 
and g\ as in (|24|) and $2 an even polynomial function of degree 12 of the form 



g2(y) = l + J2 
fc=i 



a 2 ky 2k , 



to be determined. By statement (ii) of Lemma 14. 2 \ the function M2, given in (|25p . 
associated to this V2 and k = 2/3 is of the form M2 = 4>(y)x + ip(y), where (f>(y) and 
ifj(y) are polynomials in the variable y whose coefficients depend on the unknowns a2k 
with k = 1 . . . 6. 

We fix 04 and a^ in such a way that 4>{y) has a zero at y = of multiplicity nine; we 
get the value of a$ by imposing that V2 vanishes at the two saddle points; the values 
of 02 and aio are chosen so that the curve V2 = is tangent to both separatrices at 
the saddle points of the system. Finally, after experimenting with several values for 
CJ12 and m, so that the region with boundary {V2 = 0} is positively invariant, we fix 
a 12 = -157(10m - 3)(35m + 3)/44550000. 

The region 1Z will be the bounded connected component of M 2 \ {V2 = 0} containing 
the origin, see Figured] (a). 

We need to prove that the curve {V2 = 0} n S (see Figure H] (b)) is such that the 
vector field X points in on all its points. We introduce the new parameter m = n 2 and 
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(a) (b) 

Figure 4. Positively invariant region 1Z with boundary {V 2 = 0}. 



we compute V 2 = (VT^,X) and 
Res{V 2 ,V 2 ,x) 



y*{ny 2 -lf{P n A 2 {y)fP n>m {y) 



(31) 



n 28 (120n 3 + 113n 2 - 3) 6 

where P n ,\2{y) and P n ,m{y) are polynomials of degree 12 and 36, respectively, and whose 
coefficients are polynomial functions in the variable n. 

Notice that since m € [0.5,0.547] then n € T := [0.707,0.7396]. Since the denomina- 
tor of (|31|) is positive, we only need to study its numerator. 

Using once more Lemma [3.6l and the same tools that in the previous sections we prove 
that P n ,i2{y) ■ Pn, 36 (y) is positive for all y £ (— 1/n, 1/n) and n £ T. We omit the details. 

Hence, we have proved that the numerator of Res(V2, V 2 ,x) is non-negative and it only 
vanishes on y = and y = ±n -1 ' 2 . Therefore the sets {V 2 = 0} and {V 2 = 0} only can 
intersect on {y = 0}. Indeed, the sets {V 2 = 0}nSn{y = 0} and {V 2 = 0}nSn{y = 0} 
coincide and have two points (±x(n),0) for each n 6 T. Studying the local Taylor 
expansions of V 2 (x,y) and V 2 (x,y) at these points we get that the respective curves 
V 2 (x, y) = and V 2 (x,y) = have at them a fourth order contact point and, as a 
consequence, V 2 does not change sign on {V 2 = 0} n S, as we wanted to prove. That, 
on {V 2 = 0}, the vector field X points in, is a simple verification. Hence the proof 
follows. □ 

7. Existence of polycyles 

This section is devoted to prove that the phase portrait (b) in Figure [1] can only 
appear for finitely many values of m. Notice that this phase portrait is precisely the 
only one presenting a polycycle. As we have already explained, the main difficulty is 
that we are dealing with a family that is not a SCFRVF. To see that the control of the 
existence of polycycles for general polynomial 1-parameter families can be a non easy 
task, we present a simple family for which a polycycle appears at least for two values of 
the parameter. 



Example 7.1. For m = and m = 1, the planar systems 



x 

il 



-1y + (3m - A)x + (4 



2m)x 3 + xy 2 



Pm(x,y), 
m G K. 



(4 — m)x + xy' 2 — 2mx 3 — x 5 = Q m (x, y), 
have a heteroclinic polycycle connecting the saddle points located at (±\/2 — m, 0). 



(32) 



Proof. The above family has been cooked to have explicit algebraic polycycles. Consider 
the family of algebraic curves H m (x, y) = y 2 — (x 2 + m — 2) 2 = and compute 

W m (x,y) = {VH m (x,y),{P m (x,y),Q m (x,y))). 
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Doing the resultant with respect to x of W m and H m we obtain 

Res(W m (x,y),H m (x,y),x) = m 4 (l - m) A y A R(y,m), 

where R is a polynomial of degree 4 in both variables, m and y. This implies that for 
m = and m = 1 the algebraic curve H m (x,y) = is invariant by the flow of (|32j). 
These sets coincide with the invariant manifolds of the saddle points (±\/2 — m, 0) and 
contain the corresponding heteroclinic polycycles. D 



We have simulated the phase portraits of (|32p for several values of m and it seems 
that no polycycles appear for other values of m. In any case, the example shows the 
differences between SCFRVF, for which as we have discussed in Subsection 13.11 the 
polycycle usually appears for a single value of the parameter, and families that are not 
SCFRVF. 

Let us continue our study of system ([1]). We denote by p^j = (±m -1 ' 4 , ±m -1 ' 4 ) the 
two saddle points of the system. 

Proposition 7.2. Let (0,y s (m)) be the first cut of the stable manifold of p+ with the 
Oy^ -axis. Similarly, let (0,y u (m)) be the first cut of the unstable manifold o/p~ with 
the same axis, see Figured (a). Then the function 5(m) := y s {m) — y u [m) is an analytic 
function. 




2/2 («) 

V s (m) = Jfe(j/i(m)) 




:c 



(a) (b) 

Figure 5. Definition of the maps 5(m) and y s (m) in Proposition 17.21 



Proof. This result is a consequence of the tools introduced in [21]. We only give the key 
points of that proof. 

Fix a value fh for which 8{m) is defined. Simply because the Oy + is transversal for 
the flow, the function 5 is well defined in a neighborhood of fh. It is clear that it suffices 
to prove that y s (m) is analytic at m = fh, because the y u (rn) can be studied similarly. 
To prove this fact we will write the map y s (m) as the composition of two analytic maps. 

Consider a vertical straight line £ := {{x,y) : x = fh" 1 ' 4 — e}, for e > small 
enough. Denote by (m -1 ' 4 — e,yi(m)) the first cutting point of the stable manifold of 
p^ with this line. Because C is close enough to the saddle point it can be seen that the 
local stable manifold cuts this line transversally. Moreover, the tools given in |21j prove 
that yi(m) is analytic at m = fh, because of the hyperbolicity of the saddle point. Next, 
consider the orbit starting on C with y-coordinate y\{fh). In backward time, this orbit 
cuts also transversally the Oy + -axis at the point with y-coordinate y s (fh) and needs a 
finite time to arrive to this point see Figure [5] (b). Because of the transversality to both 
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lines, and the finiteness of the needed time for going from one to the other, it is clear 
that the map 2/2(2) induced by the flow of the system between C and the Oy + -axis is 
analytic at z = y\(fh). Since y s (m) = y2(,yi(jn)), the result follows. □ 

Proof of (Hi) of Theorem \1.1[ Notice that each value of m that is a zero of the map 
5(m), introduced in Proposition 17.21 corresponds to a system (JTJ) with a polycycle, i.e. 
M = {m £ (0.547,0.6) : 6(m) = 0}. From Proposition O we know that 5(0.547) > 
and from Proposition 14.41 that 5(0.6) < 0. Hence the set M. is non-empty. Finally, 
because of the non-accumulation property of the zeros of analytic functions, the finiteness 
of M follows. □ 



8. Proof of Theorem 11.11 

The proof of Theorem 11.11 simply consists in gluing the corresponding results proved 
along the paper. More concretely: 

• The non existence of limit cycles and polycycles when m G (—00, 0.547] U [3/5, 00) 
is given in the following results: 

— For m G (—oo,0], trivially in the introduction. 

— For m G (0,9/25] in Proposition 14. 1\ 

— For m G (9/25, 1/2) in Proposition 16.11 

— For m G [1/2,0.547] in Proposition E21 

— For m G [3/5, 00) in Proposition 14.41 

• The existence of at most one limit cycle and one polycycle when m G [1/2, 3/5), 
the fact that they never coexist, and the hyperbolicity and instability of the limit 
cycle, in Proposition 15.21 

• The phase portraits of the system in the Poincare disc and the study of the 
origin, in Subsection 13.21 and Section [21 respectively. 

• The proof of the existence of the phase portrait (b) in Figure [TJ only for finitely 
many values of m, in Section [71 
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